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Abstract — Let H{D) be the parity-check matrix of an LDPC 
convolutional code corresponding to the parity-checli matrix H 
of a QC code obtained using the method of Tanner et al. We 
see that the entries in H{D) are all monomials and several rows 
(columns) have monomial factors. Let us cyclically shift the rows 
of H. Then the parity-check matrix H'{D) corresponding to the 
modified matrix H' defines another convolutional code. However, 
its free distance is lower-bounded by the minimum distance of 
the original QC code. Also, each row (column) of H'{D) has a 
factor different from the one in H{D). We show that the state- 
space complexity of the error-trellis associated with H'{D) can 
be significantly reduced by controlling the row shifts applied to 
H with the error-correction capability being preserved. 

I. Introduction 

In this paper, we assume that the underlying field is F = 
GF(2). Let G{D) be a polynomial generator matrix for an 
{no,kQ,iy) convolutional code C with memory i'. Denote by 
H(D) a corresponding parity-check matrix. Both G{D) and 
H{D) are assumed to be canonical [4], [5]. In this case, the 
code-trellis module associated with G{D) and the error-trellis 
modules [6] associated with the syndrome former H^{D) (T 
means transpose) have 2" states, where the obvious realization 
of G{D) and the adjoint-obvious realization [3] of H^{D) 
are assumed, respectively. Ariel and Snyders [1] presented 
a construction of an error-trellis based on the scalar check 
matrix derived from H{D). They showed that when some 
(jth) "column" of H{D) has a factor (i.e., the jth column 
is not "delay free"), there is a possibility that state-space 
reduction can be realized. Being motivated by their work, we 
also examined the same case. We took notice of a syndrome 
generation process. The time-fc error and the time-A: syn- 
drome are connected with the relation Ck = ^kH'^{D). 
From this relation, we noticed [8] that the transformation 



e)!_i is equivalent to dividing the jth column 



of H{D) by D . That is, reduction can be accomplished by 
shifting the "subsequence" {e^f''} of the original error-path. 



On the other hand, consider the parity-check matrix 

1 

1 1 + D + D'^ 



(1) 



Since Hi{D) is canonical and all the columns are delay free, 
any further reduction seems to be impossible. In fact, it follows 
from Theorem 1 of [1] that the dimension di of the state 



space of the error-trellis based on Hj{D) is 4. However, a 

corresponding generator matrix is given by Gi{D) = {1 + D + 
D"^, l,D^ + D^). Note that the third column of Gi{D) has 
a factor D^. {Remark: It suffices to divide the third column 
by in order to obtain a reduced code-trellis.) This fact 
implies that a reduced error-trellis can be constructed [1], [8] 
(i.e., state-space reduction can be realized). Then consider the 
reciprocal dual encoder [4] 





1 



1 



l + D + D^ 



(2) 



Note that the third column of Hi{D) has a factor D^. 
Accordingly, dividing the third column of Hi{D) by D^, we 
can construct an error-trellis with 4 states (i.e., di — 2) [1], 
[8]. Here, notice that each error-path in the error-trellis based 
on Hi{D) can be represented in time-reversed order using 
the error-trellis based on H^{D). Hence, a factor in the 
column of Hi{D) corresponds to backward-shifting by two 
time units (i.e., D~^) in terms of the original Hi{D). Actually, 
by "multiplying" the third column of Hi{D) by D^, we have 

L>2 £(2 

' (3) 



A 



1 



1 + D + D^ 







Note that this matrix can be reduced to an equivalent canonical 
parity-check matrix 



A 



1 

l + D + D^ 



(4) 



by dividing the first "row" by D^. Hence, the dimension di 
can be reduced to 2. (Remark: This fact cannot be derived 
from the results of [1].) It follows from the above argument 
that there is a possibility that a reduced error-trellis can be 
constructed not only using forward-shifted error subsequences 
but also using backward-shifted error subsequences. 

Now, we remark that a parity-check matrix H{D) with 
the form described above appears in [10]. Tanner et al. [10] 
presented a class of algebraically constructed quasi-cyclic 
(QC) LDPC codes and their convolutional counterparts. It 
is stated that the convolutional codes obtained in the paper 
typically have large constraint lengths and therefore the use 
of trellis-based decoding is not feasible. However, the parity- 
check matrices of LDPC convolutional codes proposed by 



Tanner et al. have monomial entries. Accordingly, the above- 
mentioned state-space reduction method can be directly ap- 
plied to those parity-check matrices. Then we intended to 
evaluate the state-space complexity of the error-trellis of an 
LDPC convolutional code which appears in [10]. We show 
that the overall constraint length (abbreviated as "OCL" in 
this paper) of the parity-check matrix which specifies an 
LDPC convolutional code can be significantly reduced with 
the error-correction capability of the convolutional code being 
preserved. 

II. Preliminaries 

A. Error-Trellis Construction Using Shifted Error/Syndrome 

Subsequences 

Let H{D) be a parity-check matrix for an (no, ko) convolu- 
tional code C. In this paper, we consider the error-trellis based 
on the syndrome former H^{D). In this case, the adjoint- 
obvious realization of H'^{D) is assumed unless otherwise 
specified. Denote by = {e^^\ ■ ■ , e^"'\ • • • , e^"°^) and 

Cfc iCi^^ • ■ ■ ■ ■ Cfc''' : ■ ■ ■ J Cfe'"'') the time-fc error and the time-fc 
syndrome, respectively, where r = uq — ko. Then we have the 
relation: 

Cfc = e.H^iD). (5) 

Assume that the ith row of H(D) has the form 

{D^^h[,{D) D'%^{D) ... D'^h[^^XD)), (6) 

where li > 1. Let H'{D) be the modified version of H{D) 
with the ith row being replaced by 

( K,iD) hr^{D) ... /i^„„(D) ). (7) 

Defining as ^ D^^C^^ = C,«,, we set C'^ = 

(Ci'\---,Cfe'^---,Cr)-Then we have 

C'fc = BkH'^iD). (8) 

Similarly, assume that the jth column of H{D) has the form 

{D^^h[^{D) D^^h'^^iD) ... D^^h'^^{D)f, (9) 

where Ij > 1. (Remark: H{D) is not basic [2] and then not 
canonical.) Let H'{D) be the modified version of H{D) with 
the jth column being replaced by 

( h',^{D) h',j{D) ... KjiD) (10) 

Also, let e' ^ (e^, ■ • • , e'^^^^ • ■ • , ei"")), where e',^^^ ^ 
DhfSj) = • Then we have 

Ck = eiH'^{D). (11) 

Noting these relations [8], [9], in the case where the ith 
row of H[D) has a factor by shifting the «th syndrome 
subsequence by h time units, whereas in the case where the 
jth column of H{D) has a factor D^^ , by shifting the jth 
error subsequence by Ij time units, we can construct an error 
trelUs with reduced number of states. In the following, we call 
factoring out £>' from a row of H{D) and from a column of 
H[D) "row operation" and "colurmi operation", respectively. 



t=0 t=1 t=2 t=3 t=4 1=5 




c<;'cf'=oi c?'e"=io s"s"=oi ci"e"=io s'cf=oo 



Fig. 1. Example error-treUis based on H^iD). 

B. Error-Trellis Construction Based on a Reciprocal Dual 
Encoder 

Consider the (3, 1, 2) convolutional code C2 with canonical 
parity-check matrix given by 

In this subsection, we discuss using this specific example. 
However, the argument is entirely general. Since the columns 
of H2{D) are delay free, the dimension d2 of the state space 
of the error-trellis based on the syndrome former i/J(Z)) is 
given by 2 (see Theorem 1 of [1]). Fig.l shows an error-trellis 
constructed based on _ffJ(_D) using the conventional method 
[6]. It is assumed that a transmitted code-path is terminated in 
the all-zero state at f = 4 and the corresponding received data 
is given by 2; = Zi Z2 Z3 24 Z5 = 010 Oil 000 001 000, 
where Z5 — 000 is the imaginary received data. Let z be 
the input of the syndrome former iJj(L'), then we have the 
syndrome sequence C = Ci C2 C3 C4 C5 = 01 10 01 10 00. 
The overall error-trellis is constructed by concatenating five 
error-trellis modules corresponding to Cfc- Note that the error- 
trellis in Fig.l is terminated in state (00) at t = 5, which 
corresponds to the final syndrome-former state 0-5 = (00). 
From Fig.l (note that 65 = 000), we have four admissible 
error-paths: 

= 010 Oil 000 001 000 (13) 
= 010 101 101 000 000 (14) 
= 100 000 100 000 000 (15) 
ep^ = 100 110 001 001 000. (16) 

Next, consider the reciprocal dual encoder 

Let z = Z4 Z2 Zi Zq = 001 000 Oil 010 000 be the 
time-reversed received data of {zk}\^i augmented with the 
imaginary data zq = 000. If z is inputted to the syndrome 
former ffJ(Z?), then the time-reversed syndrome sequence 
C = C5 C4 C3 C2 Ci = 00 10 01 10 01 is obtained. The 




a"ft"=oo ci"S'=io ;i,"c?'=oi (;i"??=io cS"cf'=oi 



Fig. 2. Error-trellis based on H^iD). 



corresponding error-trellis is shown in Fig. 2, where the trellis 
is terminated in state (00), which corresponds to the final 
syndrome-former state = (00). From Fig.2, we have four 
admissible error-paths: 



— 
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001 000 oil 010 000 

000 101 101 010 000 

000 100 000 100 000 

001 001 no 100 000. 



(18) 
(19) 
(20) 
(21) 



Compare these error-paths with those in Fig.l. We observe 
that each error-path in Fig.2 (restricted to the section [0, 4]) is 
represented in Fig.l in time-reversed order. That is, the original 
error-paths can be represented using the error-treUis associated 
with the corresponding reciprocal dual encoder. 

On the other hand, dividing the third column of H2 (-D) by 
D, we have the reduced canonical parity-check matrix 





l + D 



(22) 



In this case [1], [8], error-paths associated with iJj(Z)) can 
be represented using the error-trellis constructed based on 
H'2^{D). Note that a factor £)' in the column of H{D) 
corresponds to backward- shifting by / time units (i.e., D"') 
in terms of the original H{D). This observation implies that 
error-trelhs state-space reduction can be equally accomphshed 
using backward- shifted error subsequences. 

III. QC Codes and Corresponding LDPC 

CONVOLUTIONAL CODES 
A. LDPC Convolutional Codes Based on Circulant Matrices 

Each circulant in the parity-check matrix of a QC block code 
can be specified by a unique polynomial; the polynomial repre- 
sents the entries in the first column of the circulant matrix. For 
example, a circulant matrix whose first column is [1 1 1 1 0]^ 
is represented by the polynomial 1 + D + + Using this 
correspondence, an LDPC convolutional code is constructed 
based on a parity-check matrix H of a given QC code [10]. 



For example [10], let 



H = 



h 


h 


h 


h 


he 


h 


ho 


I20 


h 


/18 


I25 




I7 


hi 


hs 



(23) 



be the parity-check matrix of a (155, 64) QC code (m = 31), 
where X,. is a 31 x 31 identity matrix with rows shifted 
cyclically to the left by x positions. From H we obtain the 
following parity-check matrix with polynomial entries: 

D Jji £,8 £)16 

H{D) = I D'' £)20 ^9 ^18 

JJ25 ^19 ]JU ]j28 



(24) 



(Remark: It is stated [10] that the LDPC convolutional code 
is obtained by unwrapping the constraint graph (i.e., Tanner 
graph) of the QC code.) Note that the polynomials in H(D) 
are all monomials. In this paper, we discuss exclusively 
using this specific example. However, the argument is entirely 
general. 

B. Reordering Rows of H and the Corresponding H{D) 

Again, consider the parity-check matrix H of the (155, 64) 
QC code. Let us cyclically shift the first block of m = 31 
rows above by one position, the middle block of 31 rows by 
five positions, and the last block of 31 rows by 25 positions. 
The resulting matrix is given by 

/ 7o h h h /l5 \ 

H'=\ h h h5 h hs ■ (25) 

\ -^0 I25 hs ho h I 

Clearly, the QC block code and its associated constraint graph 
are unaffected by these row shifts. However, the convolutional 
code obtained based on the above procedure has the parity- 
check matrix 

^13 



U'(D) = 





(26) 



We see that H'{D) is not equivalent to H{D). Two con- 
volutional codes specified by H{D) and H'{D) are in fact 
different. We also remark that H{D) and H'{D) have different 
monomial entries and accordingly, when row/column factors 
are factored out, the resulting matrices have different OCLs. 

On the other hand, we have the following important fact 
[10]: 

Property: The LDPC convolutional codes obtained by un- 
wrapping the constraint graph of the QC codes have their free 
distance dfree lower-bounded by the minimum distance dmin 
of the corresponding QC code. 

It is shown that the QC code associated with H has a 
minimum distance dmin = 20. Then dj^ee of the convolutional 
code C specified by H{D) is lower-bounded by dmin = 20. 
(It is conjectured that C has a free distance dfree of 24 [10].) 
From the above property, we also have d'^^^^ > dmin = 20, 
where d'^^^^ is the free distance of the convolutional code C 
specified by H'{D). In general, let H'(D) be the parity-check 
matrix associated with H', where H' is the parity-check matrix 



obtained by applying cyclic shifts to the rows of each block 
of the original H. Above observations imply that the OCL of 
H'{D) can be controlled to some extent with its free distance 
d'f^^^ being lower-bounded by the minimum distance dmin of 
the QC code specified by H. 

IV. Reduction of Overall Constraint Length 

A. Row/Column Operations and Their Equivalent Represen- 
tation 

Again, take the parity-check matrix H given by Eq.(23). 
Here, let us cyclically shift the first block of 31 rows above 
by one position. Then the first block [Ii I2 I4 Is he] changes 
to [Iq Ii I3 /y /15]. That is, the subscript number of each 
entry decreases by 1. According to this change in H, the 
first row of H{D) changes from [D D'^ D^^] to 

[1 D _D^^]. In terms of the power of D, this change 

is expressed as [1 2 4 8 16] ^ [0 1 3 7 15]. In general, we 
observe that each entry decreases by one (modulo 31) when 
we cyclically shift the rows above by one position. Continuing 
this procedure (Remark: it is assumed that row operations have 
been done), we see that the first row of II{D) corresponds to 
one of the following five patterns in terms of the power of D: 



(27) 



Similarly, the second and third rows of II{D) are represented 
by 



r Pi 


= [0 1 3 7 15] 


P2 


= [30 2 6 14] 


< P3 


= [28 29 4 12] 


Pi 


= [24 25 27 8] 


[ p. 


= [16 17 19 23 0] 



' Qi 


= [0 5 15 4 13] 




r Ri = 


[0 25 13 20 3] 


Q2 


= [26 10 30 8] 




R2 - 


[6 19 26 9] 


< 0,3 


= [16 21 20 29] 


and < 


^3 = 


[18 12 7 21] 


Qi 


= [27 1 11 9] 




Ri = 


[11 5 24 14] 


Q5 


= [18 23 2 22 0] 




[ R,= 


[28 22 10 17 0], 



(28) 

respectively. Hence, when we apply cyclic shifts to the rows 
of each block of the original H, the resulting II{D) can be 
specified by a pattern [P,;, Qj, RkY" . For example, consider the 



pattern [Pi, Qi, P3J 
is given by 

H{D) = 



. The corresponding II{D) (cf. Eq.(24)) 



(29) 
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D 








1 




^15 




^13 


^18 


JJ12 


1 




^21 



Since row operations have been done, let us apply column 
operations. Then we have 

1 1 P>2 
1 Di5 1 1 I . (30) 

^18 ^11 I ^3 ^8 



H'{D) 



We see that this is equivalent to the transformation from 



1 3 7 15 
5* = ( 5 15 4 13 
18 12 7 21 



(31) 



30 35 40 45 50 55 60 65 70 75 80 85 90 
overall constraint length 



Fig. 3. Overall constraint length of a reduced H'(D). 



30 35 40 45 50 55 60 65 70 75 
overall constraint length {r\) 



to 



Fig. 4. Overall constraint length of a reduced H'{D). 



3 3 2 
S" = ( 4 15 
18 11 3 8 



(32) 



where subtraction is performed in each column of S in order 
that the minimum is equal to zero. From S", the OCL of the 
reduced H'{D) is obtained as 3 + 15 + 18 = 36. 

B. Reduction of Overall Constraint Length: Search Results 

As we have seen in the previous section, there are 5 x 5 x 
5 ~ 125 patterns in total. By applying column operations to 
each pattern, we examined the OCL /i of the corresponding 
reduced parity-check matrix II'{D). The result is shown in 
Fig. 3, where the horizontal axis represents the OCL ^ and 
the vertical axis represents its frequency. Observe that the 
minimum OCL /imm is 35, whereas the maximum OCL ^max 
is 83. That is, the values of n cover a wide range. Next, 
based on the argument in Section II-B, we examined the OCL 
rj using the reciprocal dual encoder II{D) associated with 
II{D). The result is shown in Fig.4. We have rjmin = 31 and 
'rjmax — 85. Note that in this example, the minimum OCL 
is further reduced using a reciprocal dual encoder. Moreover, 
after having applied row/column operations to each II{D) 
(denote hy^II'{D) the resulting matrix), we took its reciprocal 
version H'{D). Then again applying row/column operations 
to H'{D), we examined the OCL /i' of the resulting matrix. 
Using this method, we have further reduction with respect to 
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Fig. 5. Additional overall-constraint-length reduction. 



the value of fi. The resuh is shown in Fig. 5. We observe that 
the maximum reduction A/i(= /i — > 0) of 16 is realized. 
We remark that fimin = 35 is reduced to /i' = 34 using this 
method. 

C. Efficient Search Method 

Our aim is the reduction of the OCL of a parity-check matrix 
H{D). Since row operations have been done in a pattern 
[Pi, Qj, Rk]'^ , it is desirable for figures in each column to be 
close together (i.e., the difference Sj between the maximum 
and the minimum in the jth column is small). In this case, 
each figure in the column becomes small after the column 
operation, which finally leads to the reduction of the OCL. 
Hence, we search for a pattern in which every column has 
Sj < A, where A is a predetermined search parameter 

For example, set A = 20. Consider the pattern 



30 2 6 14 
Si= \ 27 1 11 9 
18 12 7 21 



(33) 



5j are given by 12, 12, 11, 7, and 12, respectively and remain 
within A = 20. Applying column operations, we have 



/ 12 2 6 5 
[ 9 1 11 
\ 12 7 12 

The OCL is given by 12 + 11 + 12 = 35. 
Similarly, consider the pattern [Pq,Q5, R^]^: 

16 17 19 23 
5*2 18 23 2 22 
28 22 10 17 



(34) 



(35) 



In this case, Sj are given by 12, 6, 17, 6, and 0, respectively. 
Applying column operations, we have 



17 6 
2 6 5 
12 5 8 



(36) 



Again, the OCL is given by 17 + 6 
these patterns have pinmi = 35. 



12 = 35. Observe that 



V. Conclusion 

In this paper, we examined the state-space complexity of the 
error-trellis of an LDPC convolutional code derived from the 
QC block code specified by a parity-check matrix H. Since the 
entries in the corresponding parity-check matrix H{D) are all 
monomials, we can construct a reduced error-trellis using the 
method of [1] or that of [8]. We noticed that when cyclic shifts 
are applied to the rows of H, the QC code remains unchanged, 
whereas the corresponding parity-check matrix H'{D), which 
defines another convolutional code, has row/column factors 
different from those in the original H{D). That is, the OCL 
of H"{D), where H"{D) is the matrix obtained by factoring 
out row/column factors in H'{D), varies depending on the row 
shifts applied to H. On the other hand, the free distance of 
the resulting convolutional code is still lower-bounded by the 
minimum distance of the original QC code. These facts imply 
that the state-space complexity of the error-trellis associated 
with H'{D) can be controlled to some extent with the error- 
correction capability being preserved and this is our basic idea. 
By applying our method to the example in [10], we have 
shown that the OCL of the parity-check matrix of an LDPC 
convolutional code can be significantly reduced compared to 
the average one. The LDPC convolutional codes proposed 
by Tanner et al. have large constraint lengths. Therefore, it 
is stated [10] that the use of trellis-based decoding is not 
feasible. We basically agree on this point. However, it has 
been shown that an error-trellis with much lower state-space 
complexity than we imagined can be constructed, which gives 
some prospect of trellis-based decoding. 
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